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摘 要: 针对 电容 层 析 成 像 (electrical capacitance tomography，ECT) 逆 问题 求解 的 病态 性 和 不 适 定性 ， 在 压缩 感知 
(compressed sensing，CS) 的 基础 上 ， 提 出 一 种 改进 FOCUSS 的 ECT 重建 算法 。 采 用 离散 余弦 变换 (DCT) 基 将 原始 图 像 
灰 度 信号 进行 稀疏 化 处 理 ， 在 使 用 正则 化 FOCUSS 算法 求解 的 过 程 中 引入 拟 牛 顿 法 逼近 求解 中 间 稀 疏 变 量 ， 以 提高 信 
号 重 构 的 准确 性 。 仿 真实 验 结果 表明 ， 同 LBP、Tikhonov 和 Landweber 和 FOCUSS 算法 相 比 ， 改 进 的 FOCUSS 算法 
能 够 有 效 区 分 物 场 中 的 不 同 介质 ， 改 善 图 像 过 度 平 滑 的 问题 ， 减 小 图 像 误 差 至 0.23， 提 高 图 像 相 关系 数 至 0.80， 具 有 
更 好 的 成 像 效果 ， 为 ECT 图 像 重建 算法 的 研究 提供 新 的 思路 
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Image reconstruction algorithm for electrical capacitance tomography based on improved 
FOCUSS algorithm 


Ma Min, Liu Ya’nan, Yang Tao, Xue Qian 
(College of Electronic Information & Automation, Civil Aviation University of China, Tianjin 300300, China) 


Abstract: On the basis of compressed sensing(CS) , this paper proposed an improved FOCUSS algorithm for the reconstruction 
of electrical capacitance tomography (ECT) , which aiming at the ill-conditionedness and ill-posedness of the inverse problem 
of ECT. The discrete cosine transform (DCT) basis made the grayscale signals of original images sparse. In the process of solving 
the inverse problem by employing the regularized FOCUSS algorithm, it introduced the quasi-Newton method to approximate 
and solve the intermediate sparse variables to improve the accuracy of signal reconstruction. Results of the simulation show that 
compared with LBP, Tikhonov, Landweber and FOCUSS algorithm, the improved FOCUSS algorithm can effectively 
distinguish the different media in the substance field, alleviate the over-smoothing effect, reduce the image error between the 
original image and the reconstructed image to 0.23, increase the correlation coefficient to 0.80, offer better image quality, and 
provide a new idea for the research on the algorithm for the reconstruction of ECT. 
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线性 和 病态 性 的 问题 , 典型 的 ECT 重建 算法 可 分 为 非 迭 代 算法 

和 从 代 算 法 : 非 迭 代 算 法 重建 速度 较 快 但 质量 不 高 ;迭代 算 济 
自 20 世纪 80 年 代 以 来 ， 电 容 层 析 成 像 技术 (ECT) 得 到 了 重建 图 像 质量 较 高 , 但 其 成 像 的 实时 性 较 差 B]。 此 外 ,典型 的 如 
迅速 的 发 展 ， 因 其 具有 无 辐射 、 非 侵入 、 结 构 简 单 、 成 本 低 、 建 算法 大 多 是 基于 工 范 数 的 优化 方法 , 通过 平滑 图 像 使 求解 稳 
响应 速度 快 等 优点 而 被 广泛 应 用 于 工业 、 食 品 和 医疗 等 领域 [1]。 定 , 导致 图 像 分 辨 能 力 较 差 , 成 像 质量 不 高 向 。 针 对 重建 算法 存 
ECT 系统 主要 分 为 三 个 模块 : 阵列 电极 传感器 模块 ， 主 要 负责 在 的 不 足 之 处 ,本文 将 基于 压缩 感知 (CS) 的 FOCUSS 算法 引入 
获取 敏感 场 内 介质 分 布 状 况 的 边界 信息 ; 数据 采集 模块 ， 主 要 ”ECT 图 像 重 建 中 。2006 年 , 由 美国 斯 坦 福 大 学 的 Donoho 及 
采集 不 同 电极 对 之 间 的 电容 值 ， 计 算 机 成 像 模块 ， 通 过 某 种 重 。 国 加 州 大 学 洛杉矶 分 校 的 Tao 等 人 提出 了 压缩 感知 理论 [5-9， 
建 算法 利用 采集 到 的 电容 值 反 演 物 场 中 介质 的 分 布 情况 [2] 。 理论 在 稀 玻 分 解 和 等 距 约 束 的 基础 上 ， 通 过 远 低 于 奈 奎 斯 特 采 
图 像 重 建 即 为 求解 ECT 逆 问 题 的 核心 ,存在 着 欠 定 性 、 非 羊 率 的 频率 对 原始 数据 进行 采样 ， 在 较 少 采样 数据 的 基础 上 准 


0 引言 


TT 


这 


曙 


SR 


瑟 


现 


H 


[Lally 


t 


三 


0 


淫 


r 


-> 


收 稿 日 期 : 2017-11-30; 修 回 日 期 : 2018-01-30 ”基金 项 目 : 国家 自然 科学 基金 资助 项 目 (61401466) ; 中 央 高 校 基金 资助 项 目 (3122013C007) ; 
家 自然 科学 基金 委员 会 与 中 国民 用 航 室 局 联合 资助 项 目 (U1733119) ; 民航 科技 资助 项 目 (20150220) 

作者 简介 : 马 敏 (1971-) ， 女 ， 安 徽 霍 印 人 ， 教 授 ， 主 要 研究 方向 为 无 损 检测 技术 及 新 型 传感器 ; 刘 亚 楠 (1993-) ， 女 《通信 作者 ) ， 河 北 沧州 人 ， 硕 
士 研 究 生 ， 主 要 研究 方向 为 电学 层 析 成 像 〈1445220996@qq.com) ; 杨涛 (1994-) ， 男 ， 湖 南 邵阳 人 ， 主 要 研究 方向 为 数据 挖 据 与 人 工 智能 ; 薛 倩 (1987-) ， 
女 ， 山 东 济 宁 人 ， 讲 师 ， 主 要 研究 方向 为 电学 层 析 成 像 . 


也 


201804.01403v1 


国 


IV 


chinaX 


录用 稿 


确 重建 出 原始 信号 。 

前 ， 压 缩 感知 重建 算法 主要 包括 三 类 : 贪 禁 方 法 、 上 四 松 
弛 方法 和 非 凸 方法 上 。 贪 焚 算 法 可 求解 范 数 最 优化 问题 , 它 
通过 逐次 迭代 选择 未 知 量 的 支撑 ， 最 终 可 逼近 原始 信号 ， 和 常见 
的 算法 有 匹配 追踪 (matching pursuit， MP)I 算 法 和 正 交 匹 配 追 
OMP)I 算 法 ; 凸 松弛 法 将 非 
凸 优化 转换 为 凸 优化 问题 来 求解 逼近 信号 , 主要 有 基 追 踩 (basic 
pursuit, BP)09 算 法 、 迭 代 收 缩 阔 值 (iterative shrinkage thresholing， 
IST)I 算法 和 梯度 投影 法 (gradient projection for sparse 
reconstruction，GPSR )02 等 ， 非 凸 方法 主要 以 欠 定 系统 局 灶 解 
FOCUSS)H31 和 稀 琉 贝 
叶 斯 学 习 (sparse bayesian learning，SBL)H41 为 主 ， 通过 改变 罚 函 
数 求解 目标 优化 问题 。 

本 文采 用 改进 FOCUSS 从 代 算 法 
上 ， 将 原始 图 像 信号 稀疏 化 ， 改 善 了 系统 本 身 的 病态 性 ， 采 用 
正则 化 FOCUSS 算法 对 ECT 逆 问 题 优化 模型 进行 求解 ， 同 时 
用 拟 牛 顿 法 通 近 中 间 稀 朴 变量 ， 克 服 了 数值 解法 中 矩阵 求 伪 道 
的 问题 使 得 求解 结果 更 加 准确 。 该 方法 在 ECT 图 像 重 建 中 取得 
了 良好 的 效果 。 


踪 (orthogonal matching pursuit， 


法 (focal underdetermined system solver， 


首先 在 压缩 感知 的 基础 
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1.3” 重 构 算 法 
对 式 (2) 直接 求解 可 转换 为 如 下 形式 : 


a=argminl|al, st.y=®Ya 


Xx=To 


五 


上 述 工 范 数 的 求解 已 被 证 明 是 一 个 NP 难 问题 ,需要 列 出 
x 中 所 有 非 零 项 位 置 的 线性 组 合 才 有 可 能 得 到 精确 解 ， 过 程 繁 
琐 且 难以 实现 。 因 此 ， 通 常 将 该 问题 转换 为 其 他 蔡 代 模型 进行 
求解 ， 常 用 的 重 构 算法 有 贪 梦 算法 、 凸 优化 算法 和 非 凸 算法 。 
本 文 利用 非 凸 算法 中 的 改进 FOCUSS 算法 进行 求解 。 


2 FOCUSS 算法 


2.1 基本 FOCUSS 算法 

FOCUSS 算法 是 Gorodnitsky 等 人 于 1997 年 提出 09， 基 本 
FOCUSS 算法 本 质 上 相当 于 一 种 加 权 最 小 范 数 最 小 二 乘法 
(weighted minimum norm least square method，MNLS)， 要 求 初 
始 信号 必须 为 稀 琉 信号 ， 采 用 琵 范 数 代替 工 范 数 可 将 y= Ax 
问题 转换 为 


X=argmin |x 
a 


及 y=Ax (4) 


求解 过 程 中 采用 线性 变换 x=Wa， 通过 过 代 将 解 的 能 量 局 


传统 的 信号 获取 其 采样 频率 必须 大 于 信号 最 高 频率 的 两 倍 ， 
即 满 足 奈 奎 斯 特 采 样 定理 。 而 压缩 感知 利用 信号 的 稀疏 性 ， 以 
远 远 小 于 传统 采样 方法 所 获 的 数据 量 便 可 完成 初始 信号 的 完整 
复 现 ， 突 破 了 奈 奎 斯 特 采样 的 限制 。 
1.1 稀 玻 表示 
假设 gp 是 RN 空间 的 一 组 向 量 基 ， 


下 =[g, 9y..9w], 则 任意 信号 量 ,ce R* 均 可 表示 为 


构成 基 拖 阵 


N 
X= 2 09 = Ya (1) 
1 


其 中 : w 是 x 在 囊 域 的 稀 疏 表示 , 若 w 的 非 零 元 个 数 比 x 少 得 
多 或 元 素 经 排序 后 成 指数 衰减 ， 则 称 x 是 稀 玻 的 或 可 压缩 的 。 
常见 的 稀疏 基 有 离散 余弦 基 、 离 散 小 波 变换 基 、 离 散 传 里 叶 基 、 
Chirplet 基 和 curvelet 基 等 。 
1.2 观测 德 阵 

CS 理论 指出 : 若 采 用 一 个 与 稀 跑 基业 不 相关 的 测量 矩阵 
四 <Rwxw ， 对 信号 量 * 进行 投影 得 到 观测 值 > ， 则 可 通过 求解 
优化 问题 由 y 精确 或 高 概率 地 重 构 出 x ， 表 达 式 如 下 : 

y= Dx= DYa O) 

其 中 : Acs _ py 为 传 感 矩 阵 ， 如 果 观 测 和 矩阵 @ 和 稀疏 基 囊 不 
相关 ， 即 在 较 高 概率 上 满足 约束 等 距 性 (restricted isometry 
property，RIP) 条 件 ， 常见 的 观测 矩阵 分 为 确定 性 矩阵 和 随机 算 
阵 。 随 机 和 矩阵 计算 复杂 程度 较 高 ， 在 硬件 上 难以 实现 ，Donoho 
等 学 者 已 经 证 明了 随机 高 斯 窍 阵 与 大 多 数 的 稀 玻 基 均 不 相干 。 
确定 性 矩阵 可 克服 上 述 缺 点 ， 但 在 压缩 感知 理论 中 不 具有 普 适 
性 。 


部 化 ， 从 而 得 到 能 量 集中 的 稀疏 解 。 基 本 的 FOCUSS 迭代 过 程 
如 下 : 


W, = Diag (x, 1) 
qe =(4Wk) > (5) 
X= Wegr 
其 中 : 权 值 矩阵 W 取 前 一 次 迭代 结果 x， 组 成 的 对 角 阵 ， 
(AW.)* 表示 和 矩阵 的 Moore-Penrose 广义 道 。 由 于 伪 逆 的 存在 ， 
往往 导致 能 量 分 散 ， 无 法 获得 精确 的 稀 疏 解 。 
2.2 正则 化 FOCUSS 算法 
针对 基本 FOCUSS 算法 存在 的 问题 , 国内 外 许多 学 者 提出 
将 工 , 范 数 引 入 FOCUSS 算法 ， 从 而 得 到 正则 化 FOCUSS 算法 


st. y=AWg 
” (6) 


gq =arg min la 


X=Wg 


其 中 : 0<P<1， 作 | 代表 p 范 数 ， 通 过 Lagrange 乘 子 法 推导 
迭代 求解 公式 如 下 : 
W. = Diag(x.1 *") 
A = AW. 
坡 H 一 1 (7) 
和 =4 (4 4+41) y 
Xi = Wi qx 


4 可 作为 噪声 信号 的 评价 指标 ， 通 常 在 噪声 标准 差 为 c ， 
信号 长 度 为 M 时 ,2=cyV2logM 可 以 有 效 去 除 噪声 干扰 。 在 本 
文中 未 直接 引入 噪声 ， 故 将 4 视 为 正则 化 参数 进行 处 理 。 


3 ”图 像 重 构 
16 电极 ECT 图 像 重 构 中 ， 归 一 化 的 模型 可 表示 为 
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A4=SG (8) 
其 中 :4 为 Mx1 维 的 电容 向 量 ( My =120 代表 独立 测量 数目 )， 
G 为 Nx1 维 的 灰 度 值 向 量 ( N =3228 代表 被 测 区 域 的 剖 分 单 
元 数目 )，$ 为 M x N 维 的 灵敏 度 年 阵 。 
于 测量 边界 数 远 远 小 于 剖 分 单元 数 ， 即 上 述 方程 组 (7) 是 
欠 定 的 ， 因 此 稳定 解 不 存在 或 不 唯一 。 
以 LBP 算法 获得 的 灰 度 和 矩阵 作为 初始 迭代 信号 G, ， 采 用 
离散 余弦 变换 (DCT) 基 对 重建 图 像 灰 度 信号 进行 稀疏 化 处 理 ， 
即 


砚 


Co = per (9) 

CS 理论 要 求 观测 矩阵 @ 和 稀 朴 基 吊 在 较 高 概率 上 满足 

RIP 条 件 ， 本 文 将 灵敏 度 矩 阵 s 的 各 行 向 量 按 高 斯 随机 顺序 重 

新 排列 得 到 新 的 随机 观测 矩阵 @' ， 使 其 具有 一 定 的 随机 性 099， 

将 电容 值 向 量 4 各行 以 相同 的 顺序 重新 排列 得 到 新 的 投影 向 
量 变换 之 后 的 优化 约束 问题 表示 如 下 : 


q=argmin|gl, 
gq 


Ss Y=AWg (10) 


以 万 范 数 代替 L 范 数 ， 采 用 正则 化 FOCUSS 算法 进行 求 
解 ， 文 献 [13] 已 给 出 FOCUSS 算法 的 收敛 性 证 明 ， 同 时 采用 
Lagrange 乘 子 法 将 其 变 为 无 约束 优化 问题 : 


1 2 
min 7(9) =lal, + Awa -Yl (11) 


目标 函数 J(g) 的 梯度 


VJ(qg)= Diaglla)| lq+ AAW)Y (AW) 


- (12) 
ViJ(g) =(AW) (AWg —Y) 
进一步 目标 函数 J(q) 的 近似 Hessian 算 阵 为 
VD = DiagllaO| 1+AAW) (AW) (13) 
采用 牛顿 法 进行 求解 ， 其 迭代 公式 如 下 
= 一 MYC -VIX®) (14) 
其 中 : y 为 搜索 步 长 ， 且 沿 下 降 方 向 进行 搜索 ，y =1 时， 是 一 


种 特殊 的 拟 牛 顿 法 ， 本 文采 用 拟 牛 顿 法 求解 稀疏 解 gw ， 将 式 
(12) (13) 代 入 式 (14) 中 进行 求解 : 


gy 三 g® {Diagl gq(i) a +ACAW)Y (AW)} 


(15) 
{ Diag[la(D q+ A(AW) (AWg -了 )} 
求解 上 式 可 得 : 
gD 蚜 二 (AW)Y (AWI (4)77 (16) 


1 


前 一 次 稀 玻 


其 中 :JII= DiasgllaG| 汪 ]= 殉 2? ， 丈 为 权 值 矩阵 


解 ga 迭代 构成 的 对 角 阵 组 成 , 令 AW = A ，4- 工 , 则 最 后 的 
m 


迁 代 公式 改写 为 


h 
等 : 基于 改进 FOCUSS 算法 的 电容 层 析 


中 
by 


gq" =[ml+A. AT AY 
GD = Wa 
在 ECT 模型 中 ,P=1 时 可 获得 较 好 的 效果 ,将 拟 牛 顿 法 通 
近 求 解 的 中 间 稀 疏 变 量 go 带 入 式 (8)， 便 可 求 得 灰 度 信号 的 


(17) 


改进 FOCUSS 和 从 代 算 法 的 求解 步骤 如 下 : 

a) 将 LBP 算法 获得 的 灰 度 矩阵 作为 初始 迭代 矩阵 ， 并 采用 
DCT 基 对 初始 信号 稀疏 化 。 

b) 构 建 观测 矩阵 使 其 与 稀 疏 基 满 足 RIP 条 件 。 

c) 采 用 正则 化 FOCUSS 算法 对 优化 模型 进行 求解 ， 同 时 用 
拟 牛顿 法 逼近 中 间 稀 疏 变量 g, ， 将 其 代入 w = Wo 中 求 得 经 压 
缩 感知 变换 后 的 稀 玻 量 。 


中) 设置 迭代 终止 条 件 | 


ae -ao| <zor ,车 不 满足 条 件 , 则 


返回 c) 继 续 迭 代 。 
6) 将 精确 解 o* 代 入 G= 平 


号 。 


E 后 的 图 像 灰 度 信 


Wt 
re 


0 得 到 习 


er 


4 ”仿真 验证 


本 文采 用 COMOSOL 5.5 和 MATLAB 2014a 对 ECT 系统 
进行 仿 


真 ， 正 问题 采用 COMSOL 自动 网 格 训 分 ， 痢 分 单元 为 
三 角 单 元 ， 以 气 液 两 相 流 为 例 构 建 模 型 ， 空 气 的 介 电 常数 设置 
为 1, 液体 的 介 电 常 数 设 置 为 3, 仿真 流 型 选择 核心 流 、 二 泡 流 、 
四 泡 流 和 十 字 型 ,其 中 泡 状 流 痢 分 为 3228 个 和 矩形 单元 , 十字 型 
| 分 为 812 个 矩形 单元 ， 分 别 采 用 LBP 算法 、Tikhonov 正则 
化 、Landweber 友 代 算法 、FOCUSS 算法 以 及 改进 FOCUSS 算 
法 对 各 流 型 进行 图 像 重建 ， 其 仿真 验证 结果 如 表 1 所 示 。 

同时 ， 选 取 图 像 误 差 IME(image error，IME) 以 及 图 像 相 关 
系数 R(correlation coefficient，CORR) 对 重建 图 像 进行 客观 评价 
L171。 图 像 相 对 误差 越 小 ， 图 像 相关 系数 越 大 ， 说 明 重 建 图 像 的 
质量 越 高 。 


六 下 


| (17) 
G 


YG -6)G.-0 
CORR= 一 所 (18) 
(56 -G) (G6.-0)’ 
1 i=1 


其 中 : G 为 设 定 的 介 电 常 数值 ，6 为 图 像 重 建 得 到 的 介 电 常 数 
值 ，G 和 分别 为 G 和 的 均值 , n 为 图 像 的 像素 数 。 各 流 型 
的 评价 指标 如 表 2 和 3 所 示 。 

表 1 可 知 ，LBP 算法 仅 在 核心 流 中 成 像 效 果 较 好 ， 而 在 
稍微 复杂 的 流 型 中 , 只 能 显示 大 概 位 置 而 无 法 精确 成 像 , 故 LBP 
算法 经 常用 于 迭代 类 算法 的 初 值 计算 。Tikhonov 正则 化 和 
Landweber 算法 在 二 泡 流 、 四 泡 流 和 十 字 型 中 成 像 效 果 较 差 ， 
易 对 物体 的 真实 位 置信 息 造成 干扰 。FOCUSS 算法 在 泡 状 流 中 
能 够 较 清 楚 地 显示 被 测 介质 的 位 置 和 形状 ， 但 易 受 周围 电极 影 
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响 ， 边 缘 伪 影 严 重 ， 在 复杂 流 


i 型 十 字 型 中 成 像 效 果 一 般 。 改 进 


FOCUSS 算法 虽然 在 核心 流 中 的 成 像 并 无 太 大 优势 ， 但 是 在 二 
泡 流 和 四 泡 流 中 效果 明显 优 于 Tikhonov 正则 化 和 Landweber 算 
法 , 在 FOCUSS 算法 的 基础 上 改善 了 边缘 伪 影 当场 域 流 型 较 
为 复杂 时 ， 也 能 够 清楚 地 显示 物体 的 位 置信 息 。 
长 1 图 像 重建 仿真 结果 
模型 核心 流 二 泡 流 四 泡 流 十 字 型 
仿真 原型 © 和 @ 十 | 
Landweber © 8 ect © 
表 2 重建 图 像 误 差 IME 
核心 流 二 泡 流 四 泡 流 。 ”十字 型 
LBP 0.3256 0.6067 0.6921 0.8034 
Tikhonov 0.2197 0.4721 0.3243 0.6002 
Landweber 0.2563 0.3068 0.4036 0.4656 
FOCUSS 0.2610 0.2954 0.3464 0.4417 
改进 FOCUSS 0.2342 0.2743 0.2872 0.4243 
表 3 府 图 像 相关 系数 CORR 
核心 流 二 泡 流 四 泡 流 “十字 型 
LBP 0.5856 0.3868 0.2516 0.2106 
Tikhonov 0.8097 0.4221 0.7843 0.3002 
Landweber 0.6882 0.7346 0.4281 0.5698 
FOCUSS 0.6721 0.7498 0.5968 0.6153 
改进 FOCUSS 0.7324 0.7663 0.8014 0.7098 
结合 表 2、3 中 的 评价 指标 可 知 ，LBP 算法 图 像 误 差 较 大 
而 图 像 相关 系数 较 小 ， 因 此 效果 最 差 。 但 是 LBP 作为 一 种 直接 
算法 并 未 进行 迭代 ， 具 有 较 快 的 成 像 速度 。Landweber 算法 在 
二 泡 流 中 图 像 相 对 误差 小 且 相 关系 数 高 ，Tikhonov 正则 化 在 四 
泡 流 中 效果 较 好 ， 而 在 较为 复杂 的 十 字 型 中 图 像 质 量 并 不 好 。 
改进 FOCUSS 算法 在 FOCUSS 算法 的 基础 上 减少 了 图 像 边 缘 


伪 影 的 影响 ， 
四 泡 流 和 


二 泡 流 的 医 


同 Tikhonov 了 


区 


术 世 明显 提高 ， cs 


区 


E 则 化 和 


Landweber 算法 相 


像 相 关系 数 均 较 
像 相关 系数 均 优 于 另外 四 


种 算法 


比 ， 其 


高 ， 在 十 字 型 中 的 成 像 效 


。 故 本 文 


马 敏 ， 等 


UD 


的 改进 FOCUSS 算法 在 ECT 成 像 9 较 高 的 成 像 质 量 。 


5 ”结束 语 


本 文 将 基于 CS 理论 的 正则 化 FOCUSS 算法 应 用 到 ECT 图 
像 重 构 , 由 于 CS 理论 可 以 在 较 少 观测 数据 的 情况 下 精确 重建 
原始 信号 ， 故 应 用 于 ECT 图 像 重 建 中 可 以 较 好 地 改善 ECT 
系统 的 欠 定 性 和 病态 性 。 采 用 正则 化 FOCUSS 算法 进行 求解 ， 
同时 将 拟 牛 顿 法 应 用 到 中 间 稀 疏 变量 的 逼近 求解 中 ， 可 进一步 
改善 图 像 质 量 。 仿 真 结 果 表 明 : 改进 FOCUSS 算法 在 ECT 
像 重建 中 具有 良好 的 效果 ， 可 为 后 续 ECT 算法 和 
考 。 
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